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In this paper we revise the main aspects of the Hamiltonian analogy: the fact that optical paths 
are completely analogous to mechanical trajectories. We follow Schrodinger's original idea and go 
beyond this analogy by changing over from the Hamilton's principal function S to the wave function 
9 . We thus travel from classical to quantum mechanics through optics. 
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I. INTRODUCTION 

According to Arnold Sommerfeld mechanics is "the backbone of mathematical physics" ; it is a highly elaborated 
theory emerging from Newton's laws which, in turn, summarize in precise form the whole body of experience. 

Classical mechanics describes with great precision the motion of macroscopic bodies, including astronomical bod- 
ies, except when velocities and masses involved are excessively large. In the domain of microscopic bodies: atoms, 
molecules, electrons, etc., the situation is quite different, for the laws of mechanics, as well as those of electromag- 
netism, do not explain satisfactorily the experimental observations, and are not appropriate for describing motion at 
a microscopic scale. During the first third of the 20th century, a new mechanics was developed -named "quantum 
mechanics" by Max Born in 1924- which governs the laws of motion for microscopic bodies but that coincides with 
classical mechanics when applied to macroscopic bodies. 

Quantum mechanics is often presented as a brand-new development based on axioms, some of them counterintuitive, 
and which seem to be hardly related to classical mechanics and other precedent theories. The axiomatic approach 
has indeed methodological advantages, but overlooks the fact that quantum theory is intimately related with classical 
mechanics, specially with its Hamiltonian formulation. Bohr's theory of electronic orbits employed Hamiltonian 
methods when it was realized the importance of separable systems in the formulation of the quantum conditions of 
Sommerfeld and Wilson in 1916; and also in the calculations of the Stark effect performed by Epstein the same year. 

The reinterpretation of the quantum laws given by Schrodinger, Heisenberg and Dirac also emanated from Hamil- 
tonian methods. The matrix nature of the canonical variables (q,p) was introduced by Heisenberg, Born and Jordan; 
while Dirac considered the conjugated variables as non-commutative operators. On the other hand, Schrodinger 
developed the operational point of view, and going farther than the analogy between optics and mechanics already 
established by Hamilton, turned the classical Hamilton- Jacobi equation into a wave equation. There is, therefore, 
a passage going from classical mechanics to quantum mechanics through optics, a path taken by Schrodinger one 
century after Hamilton. 

In 1831 William Rowan Hamilton imagined the analogy between the trajectory of material particles moving in 
potential fields and the path of light rays in media with continuously variable refractive indices. Because of its 
great mathematical beauty, the "Hamiltonian analogy" survived in textbooks of dynamics for almost a hundred 
years, but did not stimulate any practical application until 1925, when H. Busch explained the focusing effect of the 
electromagnetic field on electron beams in optical terms, which would inspire the quick development of electronic 
microscopy from 1928 onwards. Almost simultaneously, in 1926, Erwin Schrodinger went one step beyond and, using 
the ideas presented by Louis de Broglie in 1924, moved from geometric optics to wave optics of particles . In the 
present work we remember and reconstruct this passage. 

II. CLASSICAL MECHANICS 

A typical mechanical system, to which many physical systems are reduced, consists in a collection of mass points 
interacting between them following well-known laws. Experience shows that the state of the system is entirely 
determined by the assembly of positions and velocities of all its particles. The reference system chosen to describe 
positions and velocities needs not to be Cartesian: the representation of the system can be achieved by means of 
generalized coordinates, q = <?2, ■ • • , q n ), and generalized velocities, q = (<ji, q%, ■ ■ ■ , <?„)■ The minimum number, 
n, of generalized coordinates which are needed to describe completely the state of the system is called the number of 
degrees of freedom. 

The laws of mechanics are those that determine the motion of the system by providing the time evolution of the 
positions, q(i), and the velocities , q(i), of all its components. Perhaps the most general and concise way of stating 



these laws is through a variational principle known as the principle of least action or Hamilton's principle: 




where the symbol 5 stands for an infinitesimal variation of the trajectory. Here the Lagrangian, L = L(q, q, t), is a 
function of the generalized coordinates and velocities of every particle of the system. It is a general characteristic of 
the equations of classical physics that they can be derived from variational principles. Some examples are Fermat's 
principle in optics (1657) and Maupertuis'principle in mechanics (1744). Equations from elasticity, hydrodynamics 
and electrodynamics can also be derived in the same way. 

The variational principle expressed by |T]) means that when the system moves from some initial configuration at 
a given time t\ to a final configuration at some other time t%, the real trajectory is such that the action integral (|TJ) 
reaches a stationary value; it does not matter whether it is a minimum or a maximum or a saddle point. The main 
advantage of this formulation is its independence of the system of coordinates, a property that Newton's equations of 
motion do not enjoy. 

The variational principle (fTJ) leads to the Lagrange equations: 



d dL dL 

where we have used the notation 

dL f dL dL dL 
dq \dqi' dq 2 ' ' dq n 

for the gradient of L with respect to q, and similarly for dL/dq. 

Once we know L, Lagrange equations constitute a set of n second-order differential equations in time which, 
along with fixed initial conditions, determine the time evolution of the system. For a system of particles without 
electromagnetic interactions and with masses and velocities not too large, so that relativistic effects can be neglected, 
the Lagrangian is determined in such a way that the set of differential equations Q coincides with Newton's equations 
of motion: 

dp _ _dV 
dt dq ' 

where V = V(q) is the potential energy of the system. This can be achieved if 

L = T-V, 

where T is the kinetic energy. 

In the Lagrangian formulation the state of the system is governed by a set of second-order differential equations, 
each one involving second-order derivatives of the coordinates with respect to time. In many cases, specially in general 
theoretical derivations such as the present one, it is convenient to replace Lagrange equations by an equivalent set of 
twice as many equations of first-order. The most direct way of accomplishing this would be setting q = v and then 
add these additional equations to the set ([2|), treating q and v as unknowns. 

However, a much more symmetrical formulation is obtained as follows; instead of the velocities, we introduce the 
new variables: 



called momenta. Substituting this into (J2|) we have 



I- 



Therefore, we have reduced the system of n second-order equations given by ^ to a system of 2n first-order equations 
given by ([3|) and (j4]). However, the latter system is not symmetrical because velocities appear mixed with positions 
and momenta. The next step is to introduce, instead of the Lagrangian L{q 1 q, t), a new function H(q, p,i) through 
a Legendre transformation 



H = qp L, 



(5) 



where we use the notation q • p = qkPk- Using H, it is not difficult to show that Lagrange equations take the 
symmetrical form: 

.OH dH 

q= aF' p= aq"' (6) 

which is the so called canonical form of the equations of motion or Hamilton equations. Combining the definition of 
H given by ([6]) with the Lagrange equations ([2|), we also find that 

dH _ dL 
~dt ~ ~dt' 

Thus, if L does not depend explicitly of time, then dL/dt — and consequently H = constant. In this case, we say 
that the system is conservative. 

The function H(q, p, t) is called the Hamiltonian of the system. In Newtonian mechanics, where L = T — V and T 
is a quadratic and homogeneous function of the velocities, we can sec from ^ that H = T + V and the Hamiltonian 
coincides with the energy. Furthermore, if the system is conservative, then H = E = constant, which is the law of 
conservation of energy. We should note that this discussion is only valid for inertial reference frames and not for 
accelerated ones. For non-inertial systems as, for instance, a rotating reference frame, the Hamiltonian H is constant 
but does not coincide with the energy. 

III. THE HAMILTON- J ACOBI EQUATION 

For a mechanical system with n degrees of freedom, Hamilton equations are a set of 2n first-order differential 
equations which must be solved with given initial conditions. The state of the system is specified by coordinates and 
momenta instead of coordinates and velocities of the Lagrangian formulation. Therefore, knowing the state of the 
system at a given time, the solution of the differential equations will determine the state of the system at any 
subsequent time. 

Within the Hamiltonian framework, there exists an alternative way of finding the time evolution of the system. 
This alternative is based on the theory of canonical transformations (i.e., changes of variables leaving invariant the 
form of Hamilton equations). It relies on the knowledge of certain function, S^q, i), depending on the generalized 
coordinates and time, called Hamilton's principal function. Contrary to the Lagrangian and Hamiltonian approaches, 
where many ordinary differential equations have to be simultaneously solved, Hamilton's principal function satisfies 
one single equation, a non-linear partial differential equation of first order called the Hamilton- Jacobi equation: 

dS / dS \ n 

where H is the Hamiltonian in which we have replaced momenta p by the gradient of S, V<S = dS/dq. 

The complete solution of the Hamilton-Jacobi equation is any solution depending on n arbitrary and independent 
constants a. — (ai,o:2, • • • , ot n ). It can be shown that within this formalism [H, 0], the time evolution of the system is 
given by the complete solution S(q, a, t) through its derivatives: 

d d 
p = — 5(q, a,t), (3= g^S{q, a, t), (8) 

where (3 = {(3\, (3%, • • • , n ) are another set of n arbitrary constants. Note that equations (|8|) determine, in an implicit 
way, the equations of motion 

q = q(*;«,/3) p = p(t;a,p), 

which depend on 2n arbitrary constants a and (3. 

In the case of a conservative system where the Hamiltonian does not depend explicitly on time and is a constant of 
motion, _ff(q, p) = E, the time dependence of the complete solution of the Hamilton-Jacobi equation is very simple 
and reads 

S(q,a,t) = -Et + W{q_,a), (9) 



where W(q, a) a time independent function called Hamilton's characteristic function. Substituting ^ into |(7J) we 
see that this function satisfies the following differential equation, also called Hamilton-Jacobi equation, 



where the dependence on time has been removed. 

Suppose that the mechanical system consists of a single particle of mass m, moving under the action of a potential 
V(r,t), where r = (x, y, z) are the Cartesian coordinates of the particle. In such a case the Hamiltonian reads 



H=^- + V(r,t), 



and Hamilton-Jacobi equation ((7J) takes the form 



dS 1 



dt 2m 



where 



dS 



Or 



dx 



dS 
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+ V(T,t)=0, 



dSY fdS 



dy 



Moreover in the case of a conservative field (fTTj) becomes 

\VW\ 2 = 2m[E-V{v)]. 



(11) 



'§ff=iv.sf 



(12) 



A. Geometrical solution. The Hamiltonian analogy 



Let us now analyze some geometric aspects of the Hamilton principal function S which are necessary to understand 
the Hamiltonian analogy. For simplicity, we will consider that our system consists of a single particle so that S is 
given by ([9|): 



S(r,t) = W(r) - Et, 



(13) 



and W(r) satisfies (|12p. Let us remember that the geometric representation of an equation of the form /(r) = constant, 
where / is an arbitrary function, is a surface in the ordinary three-dimensional space. Accordingly, a relation of the 
form 

S(r,t) = constant, 

in which time is involved, represents a moving surface. We will detail the time evolution of such a surface when S(r,t) 
is the Hamilton principal function (|13p. In what follows, we will not write explicitly the dependence of S on the 
constants a, since we are interested on geometric aspects which do not depend on a. 

Note that at some arbitrary instant of time to, the surface S(r, to) = C, where C is a given constant, coincides with 
the surface W(r) = C+Eto (see Fig. [T]). At a later instant to + dt the surface S(r, to + dt) — C will correspond with the 
surface W(r) — C + E(to + dt). Therefore, during the interval dt the surface S = C has moved from W(r) = C + Eto 
to W(r) = C + E(to + dt). In optics, this is exactly the way wave fronts move. In other words, from a geometrical 
point of view, the surfaces 5(r, t) = constant can be considered as wave fronts moving in ordinary space. 

We will give next a short interlude in optics where this analogy becomes even clearer. Before this, we will derive 
the velocity of displacement -the "wave velocity"- of the surfaces of constant S. Let da be the orthogonal distance 
between two surfaces separated some small time interval dt (see Fig. Q}. Then, the displacement velocity is 



v = 



da 
~dt' 



(14) 



On the other hand, if we differentiate the relation S(r, t) = W(r) — Et = C, we get dW = Edt, or equivalently 



but VW ■ dr = \VW\da. Hence, 



VW ■ dr = Edt, 



E 



\VW[ 



(15) 



c + Edt 




FIG. 1: Bidimensional representation of the time evolution of the surfaces of constant S from to and to + dto (in the figure we 
assume to — 0). da is the distance between two surfaces separated a time interval dt. 

where the value of the gradient of W is calculated through the Hamilton-Jacobi equation. In the case of a particle 
moving in a conservative field |VW| is given by Equation (fT^|) and 



E 



y/2m[E - V(r)] ' 



(16) 



showing that the "wave velocity" v is not generally uniform but depends on r. 

Let us write (fT5|) in a more suggestive way We know that in terms of S, momenta are given by p = dS/dq (see 
©)■ For a conservative system, dS/dq = dW/dq = VW and 

p = VW. 



Therefore, 



E 



(17) 



On the other hand, the gradient VS = VW is a normal vector to the surface S =constant and, as p has the direction 
of the trajectory, we can conclude that the trajectories of the particle are perpendicular to the surfaces of constant S. 
As we can see, we can obtain the trajectories by drawing the normal curves to the surfaces S = constant. Therefore, 
the mechanical trajectories behave as the light rays of optics, because the latter are perpendicular to wave fronts and 
the same happens with particle trajectories which are perpendicular to the "mechanical wave fronts" S — constant. 
This is essentially the Hamiltonian analogy between geometrical optics and classical mechanics. 



IV. AN INTERLUDE IN OPTICS 



As the title of this work indicates, the passage from classical to quantum mechanics is carried out using optics as a 
guide. As we have just seen the Hamiltonian analogy, the starting point of the journey, is precisely the resemblance 
of mechanical trajectories with light rays of geometrical optics. To proceed further we have to complete the analogy 
and establish the association between mechanics and wave optics. In this section we summarize the main aspects of 
optics that are relevant for the establishment of quantum mechanics. The reader familiar with optics may skip this 
section. 

Let 4> be an electromagnetic perturbation, that is, <p is an y component of the electric field E, or the magnetic field 
B, or the potential vector A. In certain circumstances as, for instance, in the vacuum in absence of charges or in an 
electrically neutral medium, Maxwell equations show that the perturbation <j) satisfies the wave equation 



dt 2 



(18) 



where V 2 = V • V is the Laplacian operator and v is the velocity of the wave, also known as phase velocity. 



The wave equation may have countless solutions depending on the properties of the medium, the dependence of 
the wave velocity on position and time, the boundary conditions, etc. For stationary and uniform media the wave 
velocity v is a constant independent of t and r. In this case there exists a simple and special kind of solution which 
is particularly useful for the development of optics. It is the monochromatic or harmonic plane wave: 

<f>(r,t) = Aexpji u -tj +5 j , (19) 

where A, ui and S are constants, u is a unit vector defining the direction along which the wave propagates. The 
argument of the imaginary exponential, w(u ■ r/v — t) + 5, is an angle called the phase of the wave. Note that on the 
moving planes, 

u • r — vt = constant, 

where all points have the same phase, the perturbation <f>(r,t) is constant. These plane surfaces, which move perpen- 
dicularly to u with velocity v, are called wave fronts. 

One interesting feature of the solution (|19[) is that clearly shows the double periodicity in time and space. Indeed 
</>(r, t) remains unaltered when t is changed by t + T and u • r by u • r + A, where T = 2tt/u> is the period of the wave, 
v = 1/T its frequency and A = vT the wavelenght. The refractive index is defined by 

c 

n = — 

v 

where c is the speed of light in the vacuum. The wave number k is defined by 

2vr 2tt 
k = —— = n— = nko, 
A A 

where Ao = nX and kg = k/n are respectively the reduced wavelength and the wave number in vacuum. The vector 
k = fcu is called the wave vector and it has the direction of propagation of the wave, which is orthogonal to the wave 
fronts. Finally, in terms of the wave vector the monochromatic wave (|19p can be written as 

4>(r,t) = Aexp{i (k • r - ut)} , (20) 

where we have included the factor e in the amplitude A. In an equivalent way, we may write 

(f>(r,t) — Aeyop{ik (mi • r - ct)} . (21) 

In spite of its simplicity and usefulness, the monochromatic wave is an idealization rarely found in nature. Real light 
signals are aggregates of monochromatic waves, the so-called wave trains. Moreover, many media may be dispersive 
which means that the refractive index varies with the frequency. In these circumstances real light signal take the form 
of a "wave packet" (see Appendix A for a more complete discussion on this issue) . The velocity of displacement of a 
wave packet is called group velocity, v g , and in Appendix A we show that it is given by 

In a non dispersive medium frequency and wave number are related through a linear relation, lu = kv, then v g = v 
and the group velocity coincides with the phase velocity. 



A. The eikonal equation 



Inhomogeneous media, where the phase velocity and the refractive index depend on the position, constitute another 
situation often encountered in nature. Now PU]) and (|2"Tj) are no longer solutions of the wave equation (fT5|) . Neverthe- 
less, one still looks for solutions having the double periodicity in time and space showed by the monochromatic wave. 
A relatively simple way of getting this consists in changing the constant amplitude, A, and the linear term, nu • r, of 
Equation (|2ip , by arbitrary functions of the position. Then, in an inhomogeneous medium we look for solutions of 
the wave equation in the form 



</>(r, t) = A{r) exp{ik [£(r) - ct}}, 



(23) 



where A(r) and C(r) are real functions. We will see next the conditions these functions must satisfy for Equation 
(f2"31 to be a solution of the wave equation. 
Let us note first that wave fronts are now given by 

£(r) — ct = constant, (24) 

which represent arbitrary moving surfaces depending on the form of C(r). We also observe the analogy between these 
wave fronts and the surfaces S =constant of the Hamilton-Jacobi theory: 

W(r) -Et= constant. (25) 

This ratifies our previous statement obtained by geometric considerations: the surfaces of constant S are like wave 
fronts. 

The functions A(r) and C(r) defining the inhomogeneous wave (|23p are arbitrary and they need to be known. 
However, the perturbation cj>(r, t) given by (|23|) must obey the wave equation (|18p . Plugging the former into the latter 
and separating real and imaginary parts, we get 

4 7 r 2 [|V£| 2 -n 2 ]-A 2 ^=0, (26) 



V 2 £ + 2— • V£ = 0. (27) 
A 

These equations which determine both A{r) and C(r) also ensure that the expression given in (12"3")) is a solution of 
the wave equation in an inhomogeneous medium where the phase velocity and, consequently, the refractive depend 
on the position. 

Knowing A(r) and £(r) through the exact solution of (|26|) and (|27|) is rather involved, not to say impossible. 
The search for approximated solutions is therefore necessary. The help comes from the following empirical fact: the 
electromagnetic field associated to visible light has short wavelengths, between 4 x 10~ 5 cm and 7 x 10~ 5 cm. Hence, 
one can expect a good approximation to the laws of light propagation by neglecting the wavelength. In fact, this 
procedure has proved to be quite appropriate because phenomena related to non negligible wave lengths, such as 
diffraction and interference, are apparent only after high precision experiments are performed. Geometrical optics 
is a branch of optics characterized by neglecting the wavelength. Within this approximation the laws of optics are 
expressed in the language of geometry and the energy is conveyed by means of definite curves called light rays. 

We will address the problem posed by (|26|) and (|27|) using the methods of geometrical optics. We then suppose 
that the wavelength, Ao, is small. To be more precise, we assume that Ao is much smaller than the length over which 
the refractive index substantially changes, that is, 

A |Vn| < 1. 

Taking the limit Ao — > in Equation (|26|) we obtain 

\VC\ 2 =n 2 . (28) 

The function £(r) is called the eikonal (from Greek unuov = image) and (|28p is the eikonal equation, the fundamental 
equation of the geometrical optics. Once we know £(r), the surfaces 

£(r) = constant 

define wave fronts and light rays are the curves orthogonal to them. 



B. The principle of Fermat 



The orthogonality between rays and wave fronts, the so-called "ray property" is far from being trivial. In more 
than two dimensions an arbitrary continuous family of curves cannot be considered as the orthogonal trajectories 
of some family of surfaces. In mechanics the orthogonality of trajectories to the surfaces S= constant, where S is 
the Hamilton's principal function, holds only because mechanical paths are derivable from a variational principle. 
Without the principle of least action, the ray property of mechanical paths could not be established. 



In optics, the situation is completely analogous and the ray property is a consequence of a variational principle, the 
principle of Fermat (1657). Consider a ray of light propagating in some medium with refractive index n. During the 
time interval dt, and due to the property of orthogonality, the ray has moved perpendicularly between two wave fronts 
separated da (see Fig. [T]) and has, therefore, traveled the smallest distance between successive fronts. Consequently 
the time taken to cover this distance is minimum. Any other path would result in greater lengths and longer times. 
This is Fermat 's principle of the "quickest path": 

The trajectory of light satisfies the property that if a ray of light travels between two given points of the space, it 
does so in the least possible time. 

Taking into account that 

da n 
dt = — = —da, 

v c 

where n is the refractive index, we see that minimizing t means minimizing the integral 

I=t n -da, 
Jti c 

and, since c is constant, we arrive at the principle of Fermat: 

ct 2 



jidaj = 0. (29) 

The analogy between Fermat's principle ([29|) and principle of least action ([I]) is even clearer when we write the 
latter in the form called Maupertuis' principle, which is one way of writing the variational principle for a particle 
moving in a conservative field. In this case, L = T — V and from conservation of energy we write L = 2T — E, where 
E the energy of the particle. Since E is constant, the variation of the action integral J Ldt is equal to the variation of 
the integral 2 J Tdt. On the other hand, the motion of the particle is p = mu where u — dl/dt is its velocity (dl is the 
length element of the trajectory of the particle). Thus, T = p 2 /2m = (\ / 2)p(dl / dt) and the integral to be minimized 
is 2 J Tdt = J pdl. Collecting results we arrive at the principle of Maupertuis principle: 

5 ( [ pdl] = 0. (30) 



The resemblance with Fermat's principle ((2"9"|) is eloquent. The principle (j3"U|) , a special case of Hamilton's principle 
([T]), was proposed in 1744 by Maupertuis (1698-1759) who enunciated the universal hypothesis that in every event of 
nature there is some magnitude called "action" , which is minimal. 

After this digression around the world of optics, we are in a better position to go further than the Hamiltonian 
analogy and make a great leap forward toward new physics. 



V. WAVE MECHANICS 



The first step is to realize that the eikonal equation, 

I V£| 2 = n\ 

has the same form than the Hamilton- Jacobi equation, 

|VVF| 2 = 2m[E-V{v)], 

of a particle of mass m moving in a potential field V(r). Looking at both equations we see that Hamilton's characteristic 
function W plays the same role than the eikonal C, and that the expression \j2m[E — V) plays the role of the index 
of refraction. 

This comparison suggests that classical mechanics might be viewed as the geometrical optics -i.e., the limit for 
small wavelengths- of an "undulatory motion" or "wave mechanics" where rays orthogonal to wave fronts would 
correspond to mechanical paths since the latter enjoy the ray property. This also explains the analogy between 
Hamilton's and Fermat's principles. Indeed, if we compare the variational principles given by (|29p and (|30|) we see 
that the "mechanical refractive index" 



n = v / 2m(E - V) 



should be equal or proportional to the particle momentum p. This is certainly the case, because from the conservation 
of the energy 



2m ' 

we see that 

p = v/2m(£ - V), 

which completes the analogy between Hamilton's and Fermat's principles. 

Accepting that classical mechanics is the limit for small wavelengths of a hypothetical wave mechanics, two questions 
arise at once: (i) which are the wavelengths associated to the wave motion? and (ii) what is the "mechanical wave 
equation" approaching the Hamilton- Jacobi equation for vanishing wavelengths?. We will answer the first question 
in this section and postpone the second for the next. 

Let us note first that mechanics, viewed as a sort of geometrical optics, cannot account for events depending on the 
wavelength such as interferences and diffractions. We can, nonetheless, speculate about the form of a wave equation 
that would result in the Hamilton- Jacobi equation as a limit case. However, the correspondence between classical 
mechanics and wave mechanics is not biunivocal and demands conjectures, in the same manner as geometrical optics 
is contained in wave optics but not vice versa. 

Having these considerations in mind, let us try to answer the first question posed above. If Hamilton's characteristic 
function W is the analogue of the eikonal C then, recalling that Hamilton's principal function S is related to W by 

S = W-Et, 

and taking into account that the phase of an optical wave is proportional to C(r) — ct (cf. (|23p). we make, following 
Schrodinger, the fundamental hypothesis that S is analogous to the phase of the mechanical wave. Specifically, we 
assume that Hamilton's principal function is proportional to the phase of the wave: 

S = K phase, (31) 

where K is an arbitrary constant to be identified later. Moreover, looking at Equation (|23p . we write the phase of the 
wave as: 

phase = k(£ - vt) = 2tt - vtj , (32) 

where k = 2tt/X is the wave number, A the wavelength, v = v/X the frequency and v the wave velocity. Hence, 
Equation |[3"T]) yields 



W(r) - Et = 2-kK 



CM 
v ' -vt 



A 

Identifying the spatial part of this identity we obtain the eikonal of the mechanical wave 

C(r) = ^W(r). (33) 

The identification of terms involving time leads to the following relation between energy and frequency 

E = 2nKv. (34) 

This equation coincides with Planck's quantum hypothesis, E = hv, if wc recognized the arbitrary constant K, 
introduced in (|3T|) . with Planck's reduced constant 

K=^ = h. (35) 

Z7T 

Having determined the eikonal and the frequency we can easily obtain the wavelength, answering thus our first 
question. We know that A = vjv, where v is the wave velocity and v the frequency. Moreover, from Equation 
we have v = E/p, where p is the momentum of the particle and E its energy. On the other hand, by combining (|34|) 
and (|35[) we see that v = E/h. Collecting results we arrive at, 

A = -, (36) 
P 



which is the celebrated de Broglie's wavelength. 

One may argue, following the historical course of events, that the wavelength (f3"6f was proposed by Louis de Broglie 
in 1924, two years before the birth of wave mechanics occurred in 1926. Therefore, de Broglie could not have obtained 
his wavelength in the way just explained which is based on Schrodinger's assumption (|3ip . 

Indeed, the dual nature of light was conjectured by Einstein in 1905. Light propagates as an electromagnetic wave 
but it interacts with matter as if its energy were concentrated in "packets" , the quanta of the energy. This conjecture 
was experimentally verified in 1916 by Millikan's high precision measurements on the photoelectric effect. The dual 
conception of material particles is due to Louis de Broglie who argued, in 1924, that if any wavelength had to be 
associated to a particle with momentum p in an invariant relativistic way, this could only be achieved by Equation 
(|36|) . In Appendix B we give an outline of how (|36p was originally derived by de Broglie. The work of de Broglie and 
the Hamiltonian analogy are fundamental milestones in the development of quantum mechanics. 

Summarizing what we have seen up to this point we can say that a wave mechanics, which for small wavelengths 
coincide with classical mechanics, would associate to each particle of energy E and momentum p a frequency given 
by Planck's law: 

E 

V= h- 

and a wavelength given by the de Broglie's expression (|36J) . Observe that we have arrived at this conclusion only by 
assuming that Hamilton's principal function is proportional to the phase of the mechanical wave and using Planck's 
quantum hypothesis in order to identify the constant of proportionality with Planck's reduced constant h. 



VI. THE STATIONARY SCHRODINGER EQUATION 



We are now ready to answer the second question formulated at the beginning of the last section: what is the wave 
equation that for short wavelengths goes to Hamilton-Jacobi equation? 
Using (|3"Tj) and (j3"5"|) we write the phase of the mechanical wave as 

S W(r) ~ Et 

phase = — = . (37) 

n n 

Suppose, in addition, that it exists a wave field such that its intensity indicates the density of particles, in the same 
way that the intensity of the electromagnetic field denotes the density of photons. Furthermore, let us assume that 
it is a scalar field given by a single function \&(r, i) which, according to its supposed wave properties, we write in the 
form 

*(r,t) = *o(r)e lS(r <* )/r \ (38) 
Applying equation |3T|) we can separate the spatial and temporal dependence of the field and write 

tf(r,t) = iP(r)e- iEt/h , (39) 

where 

^(r) = y (r)e lW(r)/h . (40) 
Let us now impose that ^(r,t) obeys the wave equation 

H = - 2 v 2 *. ( 41 ) 

Substituting the expression for \I>(r, t) given by (j3"9"|) we see that tp(r) satisfies the equation 

E 2 

v 2 V 2 il> + -rril> = 0. (42) 
nr 

Replacing into this equation the expression (fTH|) of the velocity of the wave associated to a particle moving under a 
potential V(r), 

E 

v /2m(l! - V) ' 



we see that tjj(r) satisfies 



h 2 

—V 2 ip+{E-VU = Q. (43) 
2m 

This is the Schrodinger wave equation for a free particle in a scalar potential field. As it is a time-independent 
equation, we can interpret it as describing the stationary, e.g. periodic, motion of a particle in a field of force. But 
we could likewise apply it to the stationary beams formed by many particles, one after another and under identical 
conditions, as they appear, for instance, in electron optics. Whatever may be the case, it seems reasonable to suppose, 
according to Bom's interpretation, that the intensity of the mechanical wave, which is proportional to the square of 
the absolute value of the amplitude 

|*| 2 = 

is proportional to the particle density at the point r = (x, y, z) measured over long times; or, what is the same in 
this case, that is proportional to the probability that the particle will be observed at this spot at any instant. More 
precisely: 

\ty\ dxdydz is proportional to the probability that a particle is observed inside some element of volume dxdydz 
centered at the point r = (x, y, z). 

Let us note, however, that the notion of waves which do not transmit energy or momentum, but which determine 
probability, had become familiar to theoretical physicists from the (unsuccessful) Bohr-Karmers-Slater theory of 1924 

We will finally mention that the potential function V(r) which occurs in the Schrodinger equation can posses 
singularities at certain points or at infinity. These points (and eventually with others) will be singular points of the 
differential equation (|43[) . Since i/j(r) is to represent a wave in the stationary state, it must be free of singularities; and 
therefore Schrodinger laid down the condition that the solution of the wave equation corresponding to a stationary 
state must be one-valued, finite and free from singularities, even at the singularities of V(r). 

As we know Equation (|4"3")) possesses solutions of this kind only for certain special values of the constant E, known 
as eigenvalues. These eigenvalues will be the only values of the energy that the particle can have when it is in a 
stationary state. No wonder that the title Schrodinger gave to his first paper was precisely "Quantification as a 
problem of eigenvalues" . 



VII. GENERALIZATIONS OF THE SCHRODINGER EQUATION 

Equation (|43|) is valid for a single particle moving under a time-independent (i.e. conservative) force field. But, 
what is the Schrodinger equation for an arbitrary mechanical system, for example a system of particles, which in the 
most general case can be non-conservative? 

Let us therefore consider a non-conservative system with n degrees of freedom with a Hamiltonian of the form 

H{q, P,t) = J2 Mq, t)p k pi + V{q, t), (44) 
k,l 

where bki (k, I = 1, 2, • • • , n) is a positive definite matrix. In this case the Hamilton-Jacobi equation ([7]) reads 

dS , ,dSdS „. . n 

flt+EMq,*)^ + V(q,*) = 0. (45) 



Following Schrodinger 's hypothesis (j3Tj) , we suppose that the phase of the mechanical wave is proportional to 
Hamilton's principal function. Consequently, we assume that the wave function has the form 



where the amplitude is given by 
Differentiating (|46|) twice we get 



*(q, t) = exp {A(q, t) + iS(q, t)/h} , (46) 



If we take into account that h is very small (H = 1.05 5 x 10 34 Js) then all terms of the second member of (|47|) are 
very small compared to the last term. Hence, 



d 2 ^ _ _i_ / as as 

dq k dqi ~ K 2 \dq k dqi 
On the other hand, the time derivative of the wave function (|46[) is 



(48) 



9* 
~dt ~ 

which, considering again that h is very small, reads 



OA 
dt 



i dS 
h dt 



9* 
~dt 



i_dS 
h~dt 



(49) 



Substituting into the Hamilton- Jacobi equation (|4"5]) the expressions of (dS/dqk)(dS/dqi) and dS/dt given respectively 
by 1(15]) and (|4*9"|) . we see that *(q, t) satisfies the equation 



k.l 



1 -r r / Nt « 9* 

9(71,9(7; ft 2 ft 9t 



0, 



(50) 



which is the Schrodinger equation for this general case. 

Before proceeding further let us recall that classical mechanics and, hence, Hamilton-Jacobi equation f4"5|) are 
obtained from wave mechanics in the limit of vanishing wavelengths. But from the de Broglie's relation, A = h/p, we 
see that the limit A — > is equivalent to h — > 0. A fact that sustains the above approximations which involve the 
neglect of powers of K. One is, however, left with the impression that the foregoing derivation connecting Hamilton- 
Jacobi and Schrodinger equations is not exact, but only approximate. This conclusion is, nonetheless, erroneous: it 
can be shown that both equations are rigorously equivalent (see Ref. [ll[ for more details). 

We will next write the Schrodinger equation ([50]) in a more compact and suggestive way. We first rewrite (f50"|) as 



k.l 



<9 2 * 
dqkdqi 



V(q, *)* = ih 



9* 

at' 



or, in the equivalent way, 



5>c«.«>(?£)(?£) + "<* 



t) 



v& = ih- 



9* 
dt' 



But from the expression of the Hamiltonian given by (|44[) we can write in a formal way 



k.l 



Therefore ([5T|) is written as 



h 9 
i dq k 



^9* 



h 9 

i dqi 



+ r(q./)=//(q,^A t 



7/ 



(51) 



(52) 



(53) 



which is the general Schrodinger equation valid for any mechanical system with a general non-conservative Hamiltonian 
of the form given by (|44|) . 



One relevant remark should be made at this point since the generalization of the Schrodinger equation given by 
provides a link between wave mechanics and matrix mechanics. The latter, developed in 1925 by Heissenberg, Born, 
Jordan and Dirac, is based on abstract formulations in terms of operators. The basic relations are the commutation 
rules, which in the case of systems of one degree of freedom reduce to 



qp-pq 



(54) 



where q and p are operators corresponding to the conjugated position and momentum. 



Now, looking at the general Schrodinger equation (jS"5|) . we observe that we could have obtained that equation had 
we made the replacements 



P — >-;5-> E — >--7^ ( 55 ) 

i oq i at 

in the equation of the energy H (q, p, t) — E and thereafter operating on the wave function ^ . Note, however, that E 
is not constant and the energy equation H = E is purely formal, just a notational ease. 

We thus see that (|55|) provides the aforementioned link between wave and matrix mechanics. Schrodinger himself 
did not, at the outset of his researches, suspect any connection between his theory and matrix mechanics. He, for 
instance, wrote in 1926: "I naturally knew about his [Heisenberg's] theory, but was discouraged by what appeared 
to me as very difficult methods of transcendental algebra" [HI]. But little afterward he proved that the two theories 
were actually equivalent 0. Thus, for instance, the commutation relation (|54p becomes an obvious identity if, using 
(jST))) . we write p = (h/i)d/dq. Indeed 

After this brief excursion on matrix mechanics, we return to the generalized Schrodinger equation. For a conservative 
system ([55]) reads 

Now H = E is constant and we may look for solutions representing harmonic waves of angular frequency ui = E/h: 

*(q,t)=V(q)e- iEt / ft (57) 
Substituting into (|56|) we see that the wave amplitude ip(q) obeys the equation 

H(q,~)il, = Eil>. (58) 



i <9q 



That is, 



n 2 Y, h kM)^- + [E-V{^ = Q. (59) 

Either (|58p or (|59|) constitute the stationary Schrodinger equation for systems with n degrees of freedom. 
Finally, if the system consists of a single particle of mass m, then 

bu = hi /2m 

(Ski is the Kronecker symbol) and ([S"9"f reduces to the Schrodinger equation previously obtained (cf. (O). 



VIII. WAVE PACKETS AND CLASSICAL MOTION 



Before ending this paper we will see, through a specific example, how the wave function is related to the classical 
motion of the system. As we have seen (see also Appendix B) classical motion has to be associated not only to a 
wave but to a wave packet, i.e., a group of wave functions. It was Schrodinger who in 1926, in a paper entitled "The 
Continuous Transition from Micro- to Macro- mechanics" (reproduced in Ref. [9j), showed how to construct, for the 
harmonic oscillator, a wave packet with its amplitude very nearly concentrated at a single point which corresponds 
to the classical motion of the oscillator. 

The Schrodinger equation for the unidimcnsional harmonic oscillator of mass m and frequency to is 

lh— = - — — + -muj 2 q 2 ^. 
at 2m oq z 2 



As can be seen by direct substitution this differential equation is satisfied by 

tf(g,t) = expj-iwt/2 + C(2moj/h)e iwt q - {muj/2h)q 2 - C 2 (l - e 2iuJt ) /2V (60) 
where C is any constant. Splitting the real and imaginary parts we write 

9( a ,t)=A{q,t)e-^^ t \ (61) 

where 

A(q,t) = cxp |— mcj (q — acostot) 2 /2h\ , (62) 
a = C{2h/muo) 1 ' 2 , (63) 

and 

(p(q,t) = cut/2 + C[q(2mtu/h) 1/2 - Ccosut] sinwt. (64) 

Equation (|6ip shows that the wave function ^(q^t) represents a wave packet of amplitude A and phase ip (see 
Equation (|A3| of Appendix A). Moreover the amplitude ([62]) is a Gaussian centered at the oscillating point: 

q = a cos tot. 

Therefore, the maximum value of the amplitude of the wave packet moves as a classical oscillator of angular frequency 
to and amplitude a. The energy of such an oscillator is E = (l/2)ma 2 to 2 which, after using (|6"3"|) . reads 

E = C 2 hw. 

But, as Schrodinger had previously shown by solving the stationary equation, the energy of the quantum harmonic 
oscillator is 

E = (n+l/2)Hu 

(n = 0, 1, 2, • • • ). This sets the value of C at 

C 2 = n+l/2, (65) 

where n = 0, 1, 2, • • • is any natural number. 

On the other hand, <J"J/* is the probability density function of finding the oscillator at time t in position q. From 
([6T]) and |62]) we have 

tt* I (q - acosuit) 2 \ 

4"** = exp < -— — — — } . (66) 

n/muj i 

If h is very small (as it is) and C is very large, while a given by (|6"3"|) is finite, then the probability density, ^4'*, is 
negligibly small except when the numerator of the argument of the exponential is approximately zero; that is when 

q = a cos tot, 

which is the equation of motion of the classical oscillator. Therefore, the classical equation of motion is, with very 
high probability, the most likely state of the oscillator. Recall that we have arrived at the classical motion assuming 
a large value of C. However this implies, owing to (|65p . a large value for n. In other words, we recover the classical 
dynamics in the limit of high quantum numbers. This is precisely Bohr's correspondence principle. 



IX. SUMMARY 



The conceptual parallelism between light rays and trajectories of material particles was first noticed at the beginning 
of the 18th century by Jean Bernouilli (1667-1748), who tried to develop a mechanical theory for the refractive index 
Q. More than one century elapsed before William Rowan Hamilton (1805-1865) fully developed this idea realizing 
that the problems arising in mechanics and geometrical optics could be treated in the same unified way. 

Hamilton worked with "characteristic" and "principal" functions, both within mechanics and optics. These functions 
have the property that by simple differentiations we can find particle trajectories and optical paths. Moreover, both 
in mechanics and optics, the characteristic function satisfies the same type of differential equation, whose solution, 
with the appropriate boundary conditions, is equivalent to the solution of the equations of motion. 

Erwin Schrodinger (1887-1961) went beyond the Hamiltonian analogy. Following previous ideas from Louis de 
Broglie (1892-1987) and Albert Einstein (1879-1955), Schrodinger wondered what would be a hypothetical wave 
motion having classical mechanics as a limiting case. This would certainly imply that there is a wave associated to 
each particle. The characterization of such a wave motion is readily undertaken by assuming that the phase of the 
mechanical wave is proportional to Hamilton's principal function. Proceeding along this way, one can easily identify 
which are the wavelengths and frequencies of the mechanical wave, as well as obtaining the wave equation, which 
turns out to be the celebrated Schrodinger equation. 

Once the merging between classical mechanics and optics was completed, there appeared a new field: wave mechan- 
ics. On the other hand, from the equivalence between wave mechanics and matrix mechanics -the latter developed 
in 1925 by Werner Heisenberg (1901-1976), Max Born (1882-1970), Pasqual Jordan (1902-1980) and P. A. M. Dirac 
(1902-1984)- emerged quantum mechanics, one of the longest-range scientific revolutions of the 20th century. 
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When the refractive index of the medium is independent of the frequency, the phase velocity v is the only speed 
appearing in the problem. However, many material media are dispersive, that is to say, the index of refraction is a 
function of the frequency. In these cases v is not unique, for the phase velocity depends now on the frequency of the 
wave: different colours travel with different speeds. 

On the other hand, monochromatic waves are idealizations that cannot be completely realized in practice. Any 
light signal is a mixture of waves of different frequencies, the so-called wave trains. Even in cases where the signal is 
supposed to have been emitted by "monochromatic sources" (for instance atoms) the corresponding wave train is not 
a pure harmonic wave but a superposition of many harmonic waves of different frequencies. Therefore, in a dispersive 
medium, where each component of the signal propagates at different speeds, the concept of phase velocity loses its 
practical meaning. 

We know, however, that real light signals have finite length, in other words, they fill only finite regions of the 
space. In this situation a new concept arises: the signal velocity, as the speed of propagation of the bulk of the signal. 
Unfortunately, a light signal, as it progresses through any medium, gets distorted. A fact that hardens the way on 
how to define its velocity. 

Let us address the relatively simple case in which a given signal <j> can be represented by an aggregate of monochro- 
matic waves, all of them traveling along the same direction, and in a dispersive medium characterized by a refractive 
index, n = n(u>), depending on the frequency. Assuming without loss of generality that the common direction of 
propagation is the x axis, we write 



where k(u>) = k{—uj) and a(— u>) = a*(uj). 

By direct substitution is not difficult to see that (f>(x, t) is a solution of the wave equation. This accounts for the 
name "polychromatic waves" sometimes given to expressions like (|A1[) . On the other hand, it also implies that the 
frequency and the wave number are related by 
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APPENDIX A: LIGHT SIGNALS, WAVE PACKETS AND GROUP VELOCITY. 




(Al) 



k(uu) = n(u>)—, 



where n{ui) is the refractive index. 

A case which often appears in practical situations is that of a quasi- monochromatic wave. In this case the amplitude 
a(u>) of the aggregating waves, only differs from zero on a relatively limited range of frequencies around certain value 
wi, where the amplitude is maximum. This is the case, for instance, when a(u>) is a narrow Gaussian centered at uj\. 
Now the light signal represented by <f>(x,i) is also confined within a finite region and we speak of a "wave packet" 
(see Fig. [2]). Moreover, as far as dispersion and absorption by the medium remain small, it is possible to identify the 
velocity of the signal with the "group velocity" which we will define below (see Equation (|A5j) ). [HJ 




FIG. 2: Wave packet at a given instant of time to and one of its Fourier components. 

Suppose then a wave packet whose amplitude presents a sharp maximum at a given frequency u>\ and quickly decays 
to zero outside this maximum. If, in addition, we assume that the medium is not highly dispersive, then k(u>) will 
be a slowly varying function of the frequency uj for which we can perform the linear approximation in the integral 
appearing in (|Alj) : 

k(w) ~ hi + a(uj - u)x) H , (A2) 

here k\ = k(u)\) and a = dk/duj\i. Substituting this approximation into (|A1D and reorganizing terms we get 

(j)(x,t) ~ A(x,t)exp{i[k 1 x — wit]}, (A3) 

where 

/oo 
a{u) exp{i(cj — w\){ax — t)}du). (A4) 
-oo 

In many cases, A(x, t) is a slowly varying function. In this situations we see from (|A3[) that the group can be 
approximately considered as a wave of frequency u)\ modulated by an amplitude A which is slowly distorted. The 
traveling speed of the packet, that is, the group velocity v g , is the velocity at which the amplitude A(x,t) moves. In 
other words, v g is the speed of some observer seeing A stationary. 

Since the amplitude changes very slowly the wave group keeps approximately its shape. Moreover, we see from 
(|A4p that the dependence of A on the space x and time t is only through the combination A(ax — t). Hence, since 
A w constant, the packet moves with the speed v g = 1/a. That is (see <|A2|) ) . 



dk 



(A5) 



The relation between the phase velocity, v, and the group velocity, v g , can be obtained by differentiating the 
equation u = kv: 

, dv . , „. 

Vg = V+ dk' ( ^ 

Therefore, in a non dispersive medium where the phase velocity is independent of k, dv/dk = 0, and group and phase 
velocities coincide. 

We finally observe that the amplitude is related to the energy density carried by the wave (the square of the am- 
plitude is proportional to the energy density). Therefore, in the approximation just considered, the energy propagates 
with the group velocity because this is the speed at which the amplitude travels. Nevertheless, this is not always true. 
For example, in the so-called anomalous dispersion regions, the group velocity can overcome the speed of light or even 
become negative. In these cases the group velocity has no physical meaning. 



APPENDIX B: THE DE BROGLIE'S WAVELENGTH 



The analogy between Fermat's and Hamilton's principles aroused in de Broglie the need for a deeper understanding 
between corpuscular dynamics and wave propagation. Following this line of reasoning, De Broglie supposed in 1924 
that, with the motion of any material particle, there is associated a wave packet in such a way that the velocity of the 
particle is equal to the group velocity of the wave packet. Let m and u be respectively the mass and the velocity of 
the particle. It is also assumed that the frequency of the associated wave is given by Planck's law 

E = flu, 

where E the energy of the particle. Then, using the expression for the energy given by the special relativity, we write 

mc 2 

TVjJ 



^/l-u 2 /c 2 

Under de Broglie's assumption the velocity of the particle u equals the group velocity, (|2"2"|) . we thus have 

dio 



Hence 



U= dk- 



du) mc 2 , , 
dk = — = — - d I 1 

u hu 



-1/2 



That is, 

„ mdu 
dk = 



h{\-u 2 jc 2 fi 2 '' 

which, after recalling the rclativistic form of the momentum, 



P = 



y/1 -U 2 /C 2 ' 

can be written as 

dk=%. 
n 

Integrating we have 

p = hk (Bl) 
which, after taking into account that k = 2ir/\, leads at once to the de Broglie's wavelength (|Sfj|) : 

A = -. 

P 

Finally, the phase velocity of the de Broglie's wave is 

hE E 
v = \v = — — = — . 

p h p 

For the relativistic free particle E/p = c 2 /u, and 

_cl_ 
u ' 

which is the phase velocity of the wave associated to the particle. Note that for photons (and for other particles whose 
mass at rest is zero) v = u = c. 
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